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The deterministic Landau-Lifshitz-Gilbert equation has been used to investigate the nonlinear 
dynamics of magnetization and the specific loss power in magnetic nanoparticles with uniaxial 
anisotropy driven by a rotating magnetic field. We propose a new type of applied field, which is 
’’simultaneously rotating and alternating”, i.e. the direction of the rotating external field changes 
periodically. We show that a more efficient heat generation by magnetic nanoparticles is possible 
with this new type of applied field and we suggest its possible experimental realization in cancer 
therapy which requires the enhancement of loss energies. 
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I. INTRODUCTION 

Nonlinear dynamics of the magnetization in single¬ 
domain ferromagnetic nanoparticle systems has received 
a considerable attention due to a wide range of their ap¬ 
plicability such as magnetic resonance imaging, ultrahigh 
density magnetic data recording, spintronics, ferrofluids 
and in biomedical engineering, in particular, drug de¬ 
livery or hyperthermia [J-Q. Among many forms of 
hyperthermia, the local induction of heat via magnetic 
nanoparticles seems a fruitful strategy with promising 
preclinical results in different cancer modes [J, Q • The 
unique feature of magnetic nanoparticle hyperthermia is 
that the energy is transported in the body by means of 
an ac magnetic field. At present the clinical application 
is partly limited due to the efficacy of the heat transfer 
and the poor controllability of temperature parameters 
Q. Thus, the study of relaxation mechanisms of mag¬ 
netic nanoparticles is a very active research field, both in 
its theoretical and material-science aspects. 

Among many controllable parameters, the applied ex¬ 
ternal magnetic field is one of the most easily variable to 
increase the efficiency of heat generation. Indeed, there 
is an increasing interest in the literature to consider the 
case of rotating external magnetic field [M3 instead of 
the commonly applied alternating one, see for example 
[Til . Some of these studies compare the efficiency of the 
two types of applied fields. For example, in [13] it was 
argued that the use of an alternating-like external field is 
favourable, in particular an orthogonal synchronised bi¬ 
directional field is proposed as a most efficient heat gen¬ 
eration. Another experimental work [llj] suggests that 
the alternating and circulating applied fields produce the 
same heating efficiency in the limit of low-frequency. 

In order to clarify the above question and to look for 
the most efficient heat generation one has to take into 
account two other important effects: the role of mag¬ 
netic anisotropy which is inevitably present in magnetic 


nanoparticles and the influence of thermal fluctuations. 
For example, in ]9j, isotropic nanoparticles were consid¬ 
ered without taking into account the thermal effects and 
the alternating applied field was found to be consider¬ 
ably more favourable than the rotating one. However, 
it was shown in [l3| that the effect of thermal fluctua¬ 
tion modifies the results obtained for the alternating field 
with isotropic nanoparticles and results in a small differ¬ 
ence between the heating efficiency of the rotating and 
alternating applied fields in the limit of small frequency. 
Let us note, the inclusion of thermal fluctuations [131 ] 
only slightly modifies the findings of the rotating field 
obtained for isotropic nanoparticles without thermal ef¬ 
fects 

Magnetic anisotropy can also influence this picture, 
specially in case of the rotating applied field. An im¬ 
portant feature of the rotating external field is the pres¬ 
ence of stable steady states (precession modes). Indepen¬ 
dently of the initial positions, the magnetic moments of 
nanoparticles tend to the steady state and the dissipated 
energy can be easily calculated. However, for relatively 
large anisotropy, more than one steady states appear [3 ■ 
Thus, one can study the influence of the transition be¬ 
tween these modes on the heating HH3- For example, in 
Refs. mM the effect of strong anisotropy on the heat¬ 
ing efficiency has been studied and its enhancement is 
shown for large frequencies near the boundary of various 
regimes of forced precession (steady states) of the param¬ 
eter space in case of a rotating applied field. Although, it 
was argued in [l6| that thermal fluctuations do not mod¬ 
ify the results but the frequencies where the enhancement 
is observed are too high for a medical treatment, i.e. for 
hyperthermia. Another example where the presence of 
uniaxial anisotropy was considered is Ref. [lj] , where it 
was shown that in the low-frequency limit, an easy-axis 
anisotropy turned out to leave unchange (or slightly de¬ 
crease) the calculated loss power if no transitions between 
the various precession modes (steady states) are taken 
into account which is the case for relatively moderate 
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anisotropy. In summary, the alternating applied field is 
found to be slightly more favourable for low-frequencies 
which is suitable for hyperthermia and the small or mod¬ 
erate anisotropy does not change this picture. 

Instead of trying to decide whether an alternating or a 
rotating external magnetic field is more suitable for mag¬ 
netic nanoparticle hyperthermia, here we propose a new 
type of applied field, which is ’’simultaneously rotating 
and alternating”. In particular, we demonstrate that a 
more efficient heat generation by magnetic nanoparticles 
is possible if the direction of the rotating external field 
changes periodically and the particles exhibit a moderate 
anisotropy. We show results of this new and successful 
attempt which enhances the loss energy drastically by 
means of abrupt changes in the applied magnetic field 
(sudden change of the direction of the rotating field af¬ 
ter every circle). The change in the direction of the ro¬ 
tating applied field dislocates the magnetic moment of 
the nanoparticle out of its steady state, i.e. the preces¬ 
sion mode. Since the steady state solution corresponds 
to minimal dissipation, a more effective heating can be 
achieved if the system is out of the steady state. We argue 
that the loss energy achievable by changing the direction 
of a rotating field is worth to study as a possible tool to 
enhance the heat in hyperthermia. The calculation has 
been done for a single, simple configuration, but the stun¬ 
ning result is worth checking experimentally. Indeed, we 
propose another new type of rotating applied field which 
has a periodically alternating direction, similarly to the 
case of ” sudden change” studied in the present work but 
with a feature of being more suitable for experimental 
realisation. 

The paper is organised as follows. In Sec El we dis¬ 
cuss the deterministic Landau-Lifshitz-Gilbert equation 
in case of uniaxial anisotropy with rotating applied field 
where overall units and parameters suitable for hyper¬ 
thermia are also considered. Some known results on the 
specific loss power and loss energy obtained in the steady 
state solutions of the Landau-Lifshitz-Gilbert equation 
for oblate (positive anisotropy) particles is summaries 
briefly in Sec EH and new results for prolate (negative 
anisotropy) particles are discussed. New findings related 
to the loss energy per cycle out of the steady state are 
shown in Sec m In Sec E we propose a novel type of 
applied field which is rotating with periodically changed 
direction. Proposal for experimental realization of the 
new applied field is discussed in Sec[Vl] Finally, Sec IVIII 
stands for the summary and a detailed discussion for the 
optimised set of parameters for magnetic particle hyper¬ 
thermia is given. 


II. LANDAU-LIFSHITZ-GILBERT EQUATION 


Out of the many phenomenological equations of mo¬ 
tion for the relaxation of magnetization [19| the Gilbert 
equation [20] has proved to give the one of the most re¬ 
alistic description of the dynamics of single-domain mag¬ 


netic particles at strong damping (Ref. Q represents an¬ 
other frequently used appoach). Such a particle, being 
too small to accommodate a domain wall, can be fully 
characterized with a single vector, its magnetic moment 
m. An important feature of Larmor precession is that 
|m| = ms does not change under the influence of the 
external field, including the anisotropy field. Hence it 
is convenient to rewrite the equation of motion of the 
magnetization m of a single-domain particle in terms of 
the unit vector M = m/ms, ms being the saturation 
magnetic moment (e.g. ms ~ 10 5 A/m for single crystal 
Fe 304 ). Then the Gilbert equation reads as 


= 7o M x 
dt 


V M k + Mo 77-77 M 
dt 


(1) 


where 70 = 1.76 x 10 11 Am 2 /Js is the gyromagnetic ratio 
of the electron spin (with opposite sign), /r 0 = 47r x 1 (U 7 
Tm/A (or N/A 2 ) is the permeability of free space, V is 
the potential energy and 77 is the damping factor, both 
of them normalized for unit M. To describe the system, 
the potential energy must contain the Zeeman energy in 
the external applied magnetic field and the anisotropy 
energy [§] 


V = -moM • H ext - ^ H a M 2 z (2) 

wher M z is the z-component of the normalized magneti¬ 
zation vector and the external applied field is a rotating 
one (perpendicular to the anisotropy field) 

H oxt = H 0 (cos(wi), sin(wt), 0), (3) 

with the angular frequency oj. We define the vector H c ff, 
which contains the external applied magnetic field and 
the effect of the anisotropy of the magnetic particle 

H c ff =- VmV 7 =- {9m x , dM y , dMz)V 

Mo Mo 

= H 0 (cos(wt),sin(wt), A e ffM 2 ), (4) 

with A e ff = H a /H 0 . Clearly, if A e g > 0 ( H a > 0), the 
anisotropy will turn the magnetization towards the z- 
axis, if A c ff < 0 ( H a < 0), into the rry-plane. Giordano 
et al. [22| treated the shape anisotropy of the isotropic 
ellipsoidal particles using a different notation, where L is 
the parameter giving the deviation from spherical sym¬ 
metry. The link to our parameters is H a = (3L — l)mg. 
Although, our model does not specify the source of 
anisotropy, given the material used in hyperthermia, we 
have to deal with shape anisotropy. 

Equation 0 implies a uniaxial anisotropy with a po¬ 
tential energy — (/ro/2)L7oA e ffM 2 , which has its minima 
at \M Z \ = 1 for A e ff > 0 and M z = 0 for A e s < 0. Two 
ellipsoids of revolution have the shapes with the geome¬ 
try of this energy: for A e ff > 0 oblate (cigar-shape) and 
for A e ff < 0 prolate (lens shape). 

The Gilbert equation 0 can be rewritten in such a 
way that it has a functional form similar to the Landau- 
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Lifshitz equation (2lJ. This is called the Landau-Lifshitz- 
Gilbert (LLG) equation, 

Am = — 7 '[M x H cff ] + a'[[M x H cff ] x M], (5) 

where 7 ' = At o 7 o/(l + <a 2 ) and a ' = l'ot with the di¬ 
mensionless damping factor a = pojorjms- For example, 
a = 0.1 was chosen in 0 and a = 0.3 was used in [22| . 

Writing Eq. © in polar coordinates (M, 0 , cp) allows 
to drop the constant (M), leaving but two equations: 

ut sin 0 + ajv cos 9 cos 0 — ajvA e ff sin 9 cos 9, 

, cos 9 sin 0 

Wi COS 0 ——- + CO — CtN ——- — WiA e ff cos 9 
sin 9 sin 9 

( 6 ) 

where col = Hoj' and cxn = Hga'. Note that © is 
written in a new coordinate system rotating with the 
applied field: the azimuthal angle {ip) has been cut into 
the rotation {cot) and a measure (0) of the lag of M with 
respect to the rotation of the applied field: ip = cot — <j>. 

Knowing that in the practice [23] H 0 rts 18 kA/m we 
find that col is of the order of 10 9 Hz. In hyperthermia 
the frequency of the applied field is advised to be chosen 
between 10 5 and 5 x 10 5 Hz, so that co is four orders of 
magnitude below col- Furthermore, a at = a col- For ex¬ 
ample, a set of parameters of © typical for hyperthermia 
(with a = 0.1 and Hq = 18 kA/m), 

co = 5 x 10 5 Hz, col = 4 x 10 9 Hz, a N = 4 x 10 8 Hz, (7) 

and the dimensionless anisotropy parameter A e ff depends 
on the shape and geometry of the nanoparticle. The left 
sides of j 6 |) being derivatives of angles with respect to 
time, the units of all terms in the equations must be s _1 . 
Their dimension can be taken of introducing a dimen¬ 
sionless ’’time” t = t/to where to = 0.5 x 10 _1 °s is chosen 
in this work, thus, e.g. the dimensionless form of © 
reads 

co —y coto = 2.5 x 10 

col —t co^to = 0.2, 

cun —> cuAfto = 0.02. (8) 

Let us note, to is chosen to be in the range of the attempt 
time To used in Neel relaxation. Typical values for tq are 
between 10~ 10 and 10~ 9 seconds. 


dfl 

dt 

d0 

dt 


III. STEADY STATE SOLUTION OF THE LLG 
EQUATION 

The solution of Eq.© is shown in Fig© for a par¬ 
ticular set of (dimensionless) parameters given in the 
caption. As expected below the critical value of the 
anisotropy! A e ff | < A c , only a single attractive fixed point 
appears [h| , which corresponds to a steady state solution 


of the original LLG equation © (single periodic preces¬ 
sion mode), 

M x {t) = u x0 cos(wf) - u y0 sin(wt), 

M y {t) = u x q sin(wf) + u y o cos(wt), 

M z {t) = u z0 . (9) 

where u x o and u y o are determined by w, lol, &n and A e ff. 
The loss energy is calculated using these attractive fixed 
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FIG. 1: Orbit map in the rotating frame obtained by solv¬ 
ing the LLG equation, below the critical value of anisotropy 
(there is a single attractive fixed point in the figure). The pa¬ 
rameters are ajv = 0.1, co = —0.01, ujl = 0.2 and A e ff = 1.1. 


point solutions in the formula for the energy dissipated 
in a single cycle, 


E = p 0 m s 


dt H, 


leff 


dM 

dt 


= p 0 2TrmsH 0 {-u y0 ), 


( 10 ) 


(see also Eq.(12) in 0 ) which has the following form 
in the low-frequency co <C ajv and small anisotropy 
I A e ff | <C 1 limits 


E 


2Trp 0 msH 0 


OtjyLO 

“l + a N 


a N colco 3 

K + O 3 


(1 + 2Aeff) 


( 11 ) 


It is clear that for positive (negative) A e s the energy per 
cycle is decreased (increased) by the anisotropy albeit for 
relatively large frequencies. Since A e ff appears only at 
the next-to-leading terms, if one takes the low-frequency 
limit {co —y 0) the effect of anisotropy is irrelevant. This 
result has been discussed in 0 and can be seen in Fig© 
which shows that for |A e ff| < A c , and for small frequen¬ 
cies, the loss energy per cycle becomes identical to that 
in the isotropic case, independently of the sign of A e ff. 
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For large anisotropy, |A e ff| > A c , (if no transition be¬ 
tween the various precession modes is taken into account) 
the loss energy per cycle tends to zero for A e ff > 0 , and 
to a constant (almost identical to the isotropic value if 
uj —} 0) for A e ff < 0. Similar observation holds for the 
loss energies in the steady states which exist only above 
the critical anisotropy, see the dashed lines in Fig. [2] For 
low frequencies, the dashed lines tend to zero for A e ff > 0 
and to a constant for A e ff < 0. As to the application 
in hyperthermia, these results confirm that particles of 
lens shape provide an enhanced loss energy, but at the 
allowed frequency the effect is not significant. It is im- 
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FIG. 2: The loss energy per cycle as a function of the 
anisotropy parameter A e fr is shown for various frequencies: 
(a) uj = 0.15, (b) u> = 0.1, (c) u> = 0.05 and (d) u = 0.01. The 
calculation is based on the attractive fixed point solutions of 
© with on = 0.1, ujl = 0.2. The vertical lines indicate 
the critical value A c . Solid and dashed lines correspond to 
loss energies of steady state solutions. The upper lines (solid, 
dashed) always correspond to A e fi < 0 while the lower lines 
(solid, dashed) are related to A e fr > 0. Dotted lines represent 
the loss energy per cycle in case of rotating field where the 
direction is changed periodically. 

portant to note that frequencies (uj = 0.01 — 0.15) shown 
in Fig. [2] are too large for biomedical application. 

IV. LOSS ENERGY PER CYCLE OUT OF THE 
STEADY STATES 

Up to now, we studied the loss energy per cycle cal¬ 
culated at the steady state solution of the LLG equation 
(i.e. at the attractive fixed point of Fig. [T]). However, it 
is a relevant question to address whether one can find a 
larger loss energy in a single cycle if it is calculated out 
of the steady states. Indeed, in the 3D-plot of Fig. [3] we 
show the loss energy obtained in the first cycle of the ex¬ 
ternal field as a function of various starting points on the 
( 9 , (f >) plane. Note that in the low-frequency limit (suit¬ 
able for hyperthermia), the solution of the LLG equation 



FIG. 3: (Color online). The loss energy, E/(2Tr[iomsHo), 
obtained in the first cycle (out of the steady states) for small 
anisotropy (A e fi = 0.05) as a function of the initial conditions 
on the (6, (j>) plane. 


always tends to the attractive fixed point very rapidly, i.e. 
it reaches the fixed point with one percent accuracy in a 
quarter of a full circle. In Fig. [3] the ’’well” corresponds 
to the attractive fixed point which produces us the low¬ 
est loss energy per cycle while the ’’hill” where the loss 
energy is the maximum is related to initial conditions 
taken at the repulsive fixed point. Though, the results 
displayed in Fig. [3] are obtained numerically, an approx¬ 
imative analytic derivation is possible for the isotropic 
case (A e ff = 0). Indeed, an approximate solution to the 
LLG equation is given by Eq.(29) of [§] which reads 

M^(t) = M|q 6C cos(Ut) — M®q 6C sin(Ut) 

+ (M 50 -M| 0 pec ) e~^r, 

M v (t) = M| pec sin(flt) + M® pec cos(fii) 

+ (M^o - M s v p 0 ec ) e~^r, 

M c (t) = y/l - [M^(t)Y - [M v (t)] 2 , (12) 

where M„, M^) represents the magnetization vector 
in a particular rotated coordinate system [9j and depends 
on the initial values of the cartesian coordinates o and 
M v o linearly. (M|q 6C , M® pec depend on the parameters 

uj,ujl,o:n defined by Eq.(28) of @ and U = yj uj 2 + w 2 ). 
According to (flOl) . the loss energy per cycle has a linear 
dependence on the initial cartesian coordinates and 
M v o as well. Thus, rewriting them by polar coordinates 

o oc sin(0) cos(cj>), M v o oc sin(0) sin(</>) (13) 

one finds 

E/(lst cycle) = A sin(0) cos(^) + B sin(0) sin(</>) (14) 

where A , B constants depend on the parameters 
uj, ujl-, oln which is in agreement to the numerical result of 
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Fig. [3] obtained for small anisotropy. Finally, let us note 
that the numerical results of Fig. [3] are related to the fact 
that the steady state solution corresponds to minimal 
dissipation, so, a more effective heating can be achieved 
if the system is out of the steady state. This prelim¬ 
inary result does not define a new method to enhance 
the loss energy but the difference between the maximum 
and minimum on Fig. [3] being two orders of magnitude 
larger than what can be achieved in the steady state, is 
encouraging for further research on non-steady states. 


V. ROTATING FIELD WITH PERIODICALLY 
CHANGED DIRECTION 

As a new idea let us consider the loss energy per cycle 
in case of a rotating applied magnetic field where the di¬ 
rection of rotation is changed periodically, i.e., after every 
full cycle the sign of oj is switched in (H]). By changing 
the direction, the position of the steady state solution 
is changed, i.e. for oj > 0 (oj < 0) it is below (above) 
the equator of the orbit map. Thus, the magnetization 
vector always tends from a disappearing steady state to 
an arising one. Note that we do not calculate this ef¬ 
fect above A c because in that case more than one stable 
steady state solutions appear. The outcome is plotted in 
Fig. [2] with dotted lines. 

For high frequencies (oj = 0.15 and oj = 0.1), irrelevant 
to the sign of A e ff, there is an increase in the loss energy, 
tending to constant values close to A c . The effect is larger 
for A e ff > 0, because in this case the positions of the 
steady state solutions (for oj < 0 and w > 0) are far from 
each other. 

For low frequencies (oj = 0.001) the situation com¬ 
pletely changes as the loss energy for A e ff > 0 becomes 
larger than any other calculated values, including those 
found for A e ff < 0. For example, for lower values of the 
damping, i.e. ajv = 0.05 or ajv = 0.01 the loss en¬ 
ergy monotonously increases with A e ff, reaching a value 
that is about 15% or 100% larger as compared to the 
isotropic case near (below) A c , see Fig. [4] Therefore, a 
more efficient heat generation is possible for physical 
values of the damping (i.e. for a = 0.1 — 0.3 where 
= ctojL ~ 0.02 — 0.06), near the critical anisotropy 
A e ff = A c ~ 1 if the direction of the rotating applied 
field is changed periodically. In summary, the reason for 
the peaks of the dissipated energy at |A e s| = 1 in Fig.0] 
(dotted lines) is twofold: (i) it was shown that a more 
effective heating can be achieved if the system is out of 
the steady state which is achieved by the new type of 
applied field (dotted lines) which is simultaneously ro¬ 
tating and alternating, (ii) this enhancement effect be¬ 
comes larger at |A e ff| = 1 where the system is close to 
the critical anisotropy |A C | ~ 1, where the dislocating ef¬ 
fect caused by the change in the direction of the rotation 
is the largest. 



I \sff I 


FIG. 4: The loss energy per cycle as a function of the 
anisotropy parameter A e ff is shown for various values of qjv 
(for fixed oj and ojl)- 


VI. PROPOSAL FOR EXPERIMENTAL 
REALIZATION 

Finally, let us propose a new type of ’’rotating” ap¬ 
plied field which has a periodically alternating direction, 
similarly to the case carefully analyzed above but with 
a feature of being more suitable for experimental real¬ 
ization. Nonetheless, a ’’sharp change” in the direction 
of the rotating field which is studied previously is also 
feasible in practice. Our proposal for the applied field is 

H off = H 0 (cos(2wf),sin(wf), X eS M z ), (15) 

where the 2 oj of the cosine naturally provides us with 
the required change in the direction in every half cycle. 
Well separated and oriented anisotropic nanoparticles in 
an aerogel matrix can serve as a good experimental setup 
to realize the enhanced loss energy in case of an applied 
field as in Eq. 03 - The practical advantage of (fl5l) is 
that it requires the use of two frequencies oj and 2 oj and 
no higher harmonics or no rapid switchings are needed. 
Of course, one has to choose the frequency oj in a way 
that 2 oj should be in the range suitable for hyperthermia. 

VII. SUMMARY 

In this work we showed that a more efficient heat gen¬ 
eration by magnetic nanoparticles is possible if the direc¬ 
tion of the rotating external field changes periodically. 
We would like to emphasise that the heating efficiency 
of the new type of applied field is larger than those 
of the rotating and alternating ones which are almost 
identical in the limit of small frequency, see for exam¬ 
ple (HUH- Based on this finding, we make a proposal 
to apply a magnetic field regularly stopping its rotation, 
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which should increase the loss energy. Conditions used 
in the proposed work towards an enhance of heating ef¬ 
ficiency are the followings. 

Applied field.— A new type of ’’rotating” applied field 
is proposed here which has a periodically alternating di¬ 
rection either using the form (J4]) where w is assumed to 
change sign in every full circle or using m where the 
2 uj of the cosine naturally provides us with the required 
change in the direction in every half cycle. The frequency 
is in the range of hyperthermia, i.e. ui = 1—5 x 10 5 Hz and 
the amplitude of the applied filed is H 0 ss 0.2 x 10 5 A/m. 
Note that (fT5l) produces a change in the rotation in ev¬ 
ery half circle compared to the sudden change where the 
direction of the rotation changes in every full circle. We 
argued that in the low-frequency limit, the solution of the 
LLG equation always tends to the attractive fixed point 
very rapidly, i.e. it reaches the fixed point with one per¬ 
cent accuracy in a quarter of a full circle. Thus, it is 
expected that the proposed external field (fl5l) possesses 
the same properties found for the case of sudden change. 
Anisotropy.— The nanoparticles are assumed to be 
oblate ellipsoids where the shape anisotropy is moderate, 
i.e. A e ff = H a /Ho ss 1 (which implies L « 1.2/3 = 0.4 
according to Ref. (22| for H a = 0.2 x 10 5 A/m and 
M s = 10 5 A/m). In this case the nanoparticle has a 
single precession mode. However, for larger anisotropy, 
i.e. A c < A e ff which has been discussed in [h| (where 
h = A, , .. ), more than one precession modes appear. In 
Ref. [ly an enhancement of the heating efficiency was 
found near the boundary of these regimes of forced pre¬ 
cession. The largest increase is observed between the 
periodic and quasi periodic regimes, however, it exists 
only for large values of the reduced frequency used in [Tgj] 
which is not allowed in hyperthermia. Nevertheless, it is 
expected that in case of an applied field, Eq. (fl5l) which 
can be considered as a rotating field (with alternating di¬ 
rection) where the magnitude changes slightly, the effect 
of the increased upward heating becomes even stronger 
due to the transition between the two periodic modes. Fi¬ 
nally, let us note that regarding the anisotropy, one has 
to take into account the structure formation and syn¬ 
chronization of interacting magnetic nanoparticles which 
have been studied e.g., in [13] for rotating fields. 
Orientation.— The anisotropy field (shape anisotropy 


caused by oblate particles) is assumed to be perpendicu¬ 
lar to the plane where the applied field rotates. In prac¬ 
tice, this special situation can be achieved by switching 
on a strong static field which orients all the nanoparticles 
and than the rotating field applied while the static field 
switched off. For an appropriate size of the nanoparti¬ 
cle (with an average diameter d = 14nm [H), the ori¬ 
entation of the major part of the nanoparticles remains 
unchanged (for long enough time) and only the magnetic 
moment has a dynamics under the rotating field. Number 
of parameters (for example biomedical coating) can be 
tuned in order to design the nanoparticles (see e.g. [26j) 
for biomedical application incorporating the requirement 
for particle size and shape-anisotropy. Furthermore, in 
M it was argued that if the average diameter is chosen 
to be in the range (d m i n ,d max ) than the nanoparticles 
have single-domains and one finds no restriction on the 
applied frequency and their dynamics is almost deter¬ 
ministic. The boundary values, d m i n and d max depend 
on the parameters such as the anisotropy factor: e.g., 
in case of relatively large anisotropy (h = 0.1) one finds 
dmin = 13.7nm. Thus, our choice d = 14nm indicates 
that the description of the dynamics of the nanoparticle 
magnetic moments by the deterministic LLG equation is 
well justified [l6[ and it is expected that they retain their 
orientation for long enough time [25[. 

Damping.— The dimensionless damping constant is cho¬ 
sen to be in the range a = 0.1 — 0.3 which was used in 
the works [Til [22| . Accordingly, the enhancement effect 
shown in this work is 100% - 15% (in case of a smaller 
damping the enhancement is larger). 
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